Effects of gluon number fluctuations on 77 collisions at high energies 
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We investigate the effects of gluon number fluctuations on the total 77, 7*7* cross sections and 
the photon structure function FJ (x,Q 2 ). Considering a model which relates the dipole-dipole and 
dipole-hadron scattering amplitudes, we estimate these observables by using event-by-event and 
physical amplitudes. We demonstrate that both analyses are able to describe the LEP data, but 
predict different behaviours for the observables at high energies, with the gluon fluctuations effects 
decreasing the cross sections. We conclude that the study of 77 interactions can be useful to 



constrain the QCD dynamics. 
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I. INTRODUCTION 



The high energy evolution of dipole scattering amplitudes in QCD is described by the pomeron loop equations 
[lj-lij], a generalization of the Balitsky-JIMWLK hierarchy (see Q and references therein) by including the gluon 
number fluctuations. In the case when the strong coupling constant a e is fixed, the fluctuations wash out the BFKL 
approximation and the geometric scaling behavior predicted by Balitsky-Kovchegov (BK) equation 0-Q-the simplest 
(mean-field) nonlinear evolution equation which describes the evolution of the amplitude for the scattering between 
Q_[ a quark- antiquark pair projectile (a dipole) and a dense target. 

JL . In the last decade, many efforts have been made in the investigation of the fluctuation effects. On the theoretical 
side, because of the complexity of the pomeron loop equations, whose properties up to now have been obtained only 
under some approximations, in the last few years fluctuations in high energy evolution have been studied through 
simple (toy) models inspired in QCD |9l-[l6j. One of such models, which allows the inclusion of both fluctuation 
, and running coupling effects simultaneously, has shown that fluctuations are strongly suppressed by the running of 



the coupling, up to extremely high energies 17]. However, it should be pointed out that this result can be model 
Q\ [ dependent and that the investigation of both effects in real QCD is still a challenge. On the phenomenological side, 
f**«> ■ and in the fixed coupling case, fluctuation effects have been investigated at HERA and RHIC/LHC energies. They 
have been included in the description of the data on inclusive and diffractivc electron-proton deep inelastic scattering 
(DIS) 18-22]. Although the results have shown some improvement in the description of the observables, they have 
1 not been conclusive with respect to the presence of the fluctuations in the experimental data. They have also been 
O studied in the analysis of the pseudo-rapidity distribution of hadron multiplicities of high energy Au+Au collisions 
at RHIC and in predictions for these observables in Pb+Pb collisions by using Color Glass Condensate dynamics 
at LHC/ ALICE [23]. It has been found that the charged hadron multiplicities at central rapidity are significantly 
^ ■ smaller than saturation based calculations and are compatible to those obtained on a study of multiplicities in the 
fragmentation region with running coupling corrections [241 ] . 

One can see that, up to now, high energy QCD phenomenology in the presence of fluctuations has been studied in 
few papers only and, particularly, in processes where two scales are present. As is well known, for ep/pp colliders, the 
study of the QCD Pomeron is made difficult by the fact that the cross section is influenced by both short and long 
distance physics. Only when specific conditions are satisfied is that one can expect to determine the QCD pomeron 
effects. Some examples are the forward jet production in deeply inelastic events at low values of the Bjorken variable 
x in lepton-hadron scattering and jet production at large rapidity separations in hadron-hadron collisions, which are 
characterized by one hard scale. This motivates us to look for their effects on different processes. One of these 
processes is the off-shell photon scattering at high energy in e + e~ colliders, where the photons are produced from 
the lepton beams by bremsstrahlung (For a review see, e.g., Rcf. 25]). In these two-photon reactions, the photon 
virtualities can be made large enough to ensure the applicability of perturbative methods. Moreover, the photon 
virtualities can be varied to test the transition between the soft and hard regimes of the QCD dynamics and it is 
possible to scan the kinematical region to determine the range where the contribution of the gluon fluctuations effects 
is larger. Up to now, studies of nonlinear QCD effects on these reactions were done without taking into account gluon 
number fluctuations (see [26], [2^] and references therein) . Such investigation is the aim of the present work. 

In this paper we investigate the consequences of the inclusion of gluon number fluctuations in 7W7C*) collisions. 
Photon-photon interactions can be understood as a dilute-dilute scattering and thus favours one to rely on the 
pomeron loop equations. Within the dipole picture, we demonstrate that in the kinematical range of the LEP data on 
total 77, 7*7* cross sections and the real photon structure function F^{x, Q 2 ), the gluon fluctuations effects are small. 
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However, they contribute significantly in the range which will be probed in the future linear colliders. The paper is 
organized as follows. In Section [TT] we review some important properties of nonlinear high energy QCD evolution of 
the dipole-hadron scattering amplitude, in order to explain how fluctuation effects can be included in this analysis. 
In Section IlIII we describe the dipole representation of 77 scattering and present a model for the dipole-dipole cross 
section, the main input of the calculation of observables. Section ITVl is devoted to the description of the available LEP 
data on two photon collisions as well as predictions for future experiments. The conclusions are given in Section [V] 



II. QCD DYNAMICS AT HIGH ENERGIES 



Let us consider the general problem of a scattering between a small dipole (a colorless quark- antiquark pair) and a 
dense hadron target, at a given rapidity interval Y. The dipole has transverse size given by the vector r = x — y, where 
x and y are the transverse vectors for the quark and antiquark, respectively, and impact parameter b = (x + y)/2. 
An important quantity in the description of this process in the high energy regime is (the imaginary part of) the 
forward scattering amplitude (T(r,b)) Y = (T(x,y)) Y = (T xy ) Y , whose evolution with rapidity is given by 

cV (T xy ) Y =a J d 2 z ( x } X z y^!_ y y [( T ^)y + ( T zy) Y ~ ( T *y) Y ~ ( T v* T zy) Y ] > W 

where a = a s N c /ir and (• • •) denotes the average over all the configurations of the target at a given rapidity interval 
Y. This equation is the first equation of a infinite hierarchy, the Balitsky-JIMWLK hierarchy and has a simple 
interpretation in terms of the projectile dipole evolution: if the rapidity is increased by an ammount 5Y, there is a 
probability for a gluon, with transverse coordinate z, to be emitted by the quark (or antiquark) of the pair. In the 
large N c limit (N c is the number of colors), this gluon can be replaced by a quark- antiquark pair at point z. This 
is the dipole picture introduced by Mueller [28| . Thus, after one step in the evolution, the incoming dipole (x, y) 
splits into two new dipoles (x, z) and (z,y), which then interact with the target. The last term, (T (x , z)T '(z , y)) , 
corresponds to the scattering of both new dipoles with the target. 

If one performs a mean field approximation, (T(x, z)T(z, y)) Y ~ (T(x,z)) (T(z,y)) Y and the resulting equa- 
tion is the BK (Balitsky-Kovchegov) equation p^g], a closed equation for the average dipole scattering amplitude 
(T(x, z)) Y = A/y(r), which, at fixed coupling is given by 



d Y M Y (x,y) = a I d 2 z _ (X w y) _ , 2 Wy(x,z) + M Y {z,y) - M Y {x,y) -Af Y (x,z)Af Y (z,y)\ , (2) 
J [x z) (z y) 

This equation includes unitarity corrections and is free from the problem of diffusion to the infrared (nonperturbative) 
region (present in the solution of the linear BFKL equation [29|]). Its solution has the following properties: (i) for 
small r = \r\, Af(r) is small -the color transparency regime- and is well approximated by the BFKL solution; (ii) 
for large r, the amplitude approaches the unitarity bound Af(r) = 1, the so called 'black disc' limit, and the transition 
between these two regimes takes place at r = 1/Q S (Y). Q S (Y) is an increasing function of rapidity Y and is called 
the saturation scale, defined in such a way that Af(r) = 0(1) (usually 1/2) when r = 1/Q S (Y). 

BK equation has been shown 30] to belong to the same universality class of the Fisher and Kolmogorov-Petrovsky- 
Piscounov (FKPP) equation [31 1. Therefore, BK equation admits traveling wave solutions: at asymptotic rapidities, 
the scattering amplitude depends only on the ratio r 2 Q 2 (Y) instead of depending separately on r and Y. This scaling 
propert y is called geometric scaling and has been observed in the measurements of the proton structure function at 
HERA [32j . The amplitude is a wavefront which interpolates between and 1 and travels towards smaller values of r 2 
with speed A - the saturation exponent - keeping its shape, and the saturation scale Q S (Y) gives the front position. 

Within the correspondence between reaction-diffusion processes and the QCD evolution at high energy, it has 
been realized that the Balitsky-JIMWLK hierarchy is not complete because they do not take into account the gluon 
(dipoles) number fluctuations, which are related to discreteness in the evolution, and thus they are completely missed 
by BK equation. At least at fixed coupling, fluctuations influence dramatically the QCD evolution at high energies, 
and so the properties of the scatterin g am plitudes. Their inclusion results in a new hierarchy of evolution equations, 
the pomeron loop equations [H-U l33l. l34j . Because of the complexity of the equations of this new hierarchy, many of 
their properties have been known from some approximations [l| , after which it has been found that the hierarchy can 
be generated from a Langevin equation for the event- by- event amplitude. Formally, this is the BK equation with a 
noise term, which lies in the same universality class of the stochastic FKPP equation (sFKPP): each realization of the 
noise means a single realization of the target in the evolution and leads to an amplitude for a single event. Different 
realizations of the target lead to a dispersion of the solutions, and then in the saturation momentum p s = ln(Qg/fcg) 
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FIG. 1: Comparison between the average amplitude (Mgbw) and the event-by-event amplitude predicted by the GBW, IIM-S 
and rcBK models. The behaviours with r 2 are shown for x = 10 -3 (left panel) and 10 -6 (right panel). 



from one event to another. The saturation scale is now a random variable whose average value is given by 

(Q2(r ))= exp[A*y] (3) 

and the dispersion in the position of the individual fronts is given by 

a 2 = {pi) (p s ) 2 = DaY. (4) 

where D is the diffusion coefficient, a number expected to be of order one, which determines the rapidity Yjj = l/D 
above which gluon number fluctuations become important. 

The probability distribution of p s is, to a good approximation, a Gaussian [35[ 



Py(p s ) ^ ^=exp 



(Ps-(Ps)) 2 



(•5) 



For each single event, the evolved amplitude shows a traveling- wave pattern, which means that geometric scaling is 
preserved for each realization of the noise. However, the speed A* of the wave is smaller than the speed predicted by 
BK equation. The average (or physical) amplitude is determined by (p = ln(l/r 2 QQ)) 

/+oo 
dp s P Y {p s )M{p lPs ) , (6) 
-oo 

with Af(p, Ps ) being now the event-by-event scattering amplitude. A crucial property of the physical amplitudes is 
that at sufficiently high energies, unlike the individual fronts, they will generally not show geometric scaling. More 
specifically, they will show additional dependencies upon Y, th roug h the front dispersion a. Then, geometric scaling 
is washed out and replaced by the so-called diffusive scaling 0, 133TI361 [37| 

{ M( p ,p s )) =N (^d). (7) 

The different scaling behaviours arising from different versions of QCD evolution were compared with the available 
HERA data for inclusive, exclusive and diffractive observables in Ref. [38[ and the quality factor, which estimates 
the validity of the scaling, was determined. They found that the diffusive scaling leads to best quality factor for 
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vector meson production at HERA. Furthermore, in |18| the authors found that the description of the DIS data is 
improved once gluon number fluctuations are included and that the values of the saturation exponent and the diffusion 
coefficient turn out reasonable and agree with values obtained from numerical simulations of toy models which take 
into account fluctuations. For instance, for the event-by-event amplitude given by the GBW model [3!| 

AA GBW (r,y) = l- e -^W/ 4 , (8) 

where the saturation scale is given by Q 2 S {Y = \u(xq/x)) = Q\ (xo/x) X , they have found that A = 0.225 and D = 0.397 
for a % 2 /d.o.f. = f.f4. In contrast, for the D = case (no fluctuations), A = 0.225 and x 2 /d.o.f. = 1.74. A similar 
conclusion was obtained considering the IIM model [13] for the event-by-event amplitude, with the values of A and D 
being quite model independent. Fig. [T] shows a comparison between the behaviours for the scattering amplitude with 
r 2 for given values of x, with and without the inclusion of fluctuations. Besides GBW model for the event- by- event 
scattering amplitude we also present the predictions of the IIM-S [4l[ and rcBK [42[ models used in Ref. [13] ■ One can 
see that, when the gluon fluctuations effects are included, the onset of saturation is strongly delayed in comparison to 
the event-by-event scattering amplitude of the GBW model. The same conclusion is valid when (Ngbw) is compared 
to the IIM-S prediction. In comparison to the rcBK one, the GBW averaged amplitude has a similar r 2 behaviour 
for x — I0~ 3 . However, they become quite different at smaller values of x. 

Although a definitive conclusion is not possible from the phenomenological studies presented in Refs. [l8l - |22l [38| . 
they indicate that the presence of the gluon fluctuation effects cannot be disregarded at HERA. Moreover, the 
contribution of these effects is expected to increase with the energy, which implies that it can be large at LHC. 
However, due the complexity of the hadron-hadron collisions, it is not clear if the discrimination of these effects will 
be feasible. Consequently, the search for alternative processes to constrain the presence and magnitude of the gluon 
fluctuations effects is justified. This will be done in the following through two-photon scattering processes. 



III. PHOTON - PHOTON COLLISIONS IN DIPOLE REPRESENTATION 

Cross sections of 'yM'yM scattering can be measured at e + e~ colliders by tagging both outgoing leptons close to 
the forward direction (for a review see [25]). The process is described by the reaction e + + e~ — > e + + e~ + X, 
where A is a generic hadronic state formed through the interaction of photons emitted by the two leptons, and can 
be represented by the Feynman diagram in FigJ^J where q\ and qi are the photon four-momenta. 




FIG. 2: Feynman diagram for the e + e — > e + e X process. 



At high energies, the scattering between the two photons can be described in the dipole frame, in which the photons, 
with virtualities Q\ 2 = — q\ 2 , fluctuate into quark- ant iquark pairs (two dipoles) with transverse sizes T\ i, which then 
interact and produce the final state (see Figj3]). Within such formalism, the part of the two-photon total cross section 
that determines the energy behaviour at high energies corresponds to the exchange of gluonic degrees of freedom and 
is given by [43| 
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In the above formula, W 2 = (qi + q 2 ) 2 is the collision center of mass squared energy, are the longitudinal 
momentum fractions of the quarks in the photons, \&?(2fc,r) denotes the photon wave function, the indices i, j 
label the polarisation states of the virtual photons (i, j =L or T) and a, b label the quark flavours. a em is the 
electromagnetic coupling constant. The interaction is described by <J dd b (ri, r 2 , Y), which is the dipole-dipole cross 
section. In the eikonal approximation, it can be expressed by 



a dd {r 1 ,r 2) Y) = 2 J d 2 bAA(n,r 2 ,b,y) 



(10) 



where J\f(ri,r 2 ,b,Y) is the imaginary part of the scattering amplitude for two dipoles with transverse sizes r\ and 
r%, relative impact parameter b and rapidity separation Y . The inclusion of the unitarity corrections in the dipole 
- dipole scattering was addressed in Ref. [44| considering independent multiple scatterings between the dipole, with 
unitarization obtained in a symmetric frame, like the center-of-mass frame. Such corrections were also estimated 
considering the Color Glass Condensate formalism in Ref. 45]. As in general the applications of the CGC formalism 
to scattering problems require an asymmetric frame, in which the projectile has a simple structure and the evolution 
occurs in the target wavefunction, the use of the solution of the BK equation in the calculation of the dipole-dipole 
scattering cross section is not so straightforward. 




FIG. 3: Two-photon interactions in the dipole representation. 

In order to relate Af(ri, r 2 , b, Y) with Af(r, Y) we follow Ref. [27j |. which considers the model proposed by Iancu, 
Kugeratski and Triantafyllopoulos to study the Mueller-Navelet process [4f|. In this model (denoted IKT model 
hereafter) the dipole-dipole cross section has the following form: 

a dd (r 1 ,r 2 ,Y)=27rr 2 M(r 2 ,Y 2 )Q(r 1 -r 2 ) + 2nr 2 M(r u Y 1 )e(r 2 -r 1 ) , (11) 

where Yi = ln(l/xi) and 

* " i^Tof (12) 

The main assumptions in this model are the following: (i) The radial expansion of the gluon distribution in the target 
(larger dipole) only affects the subleading energy dependence of a dd , which implies that it is possible to study the 
approach towards unitarity limit at a fixed value of the target size; (ii) Only the range b < R, where R = Max(ri, r 2 ), 
contributes for the dipole-dipole cross section, i.e. it is assumed that AT is negligibly small when the dipoles have no 
overlap with each other (b > R). As shown in Ref. [27|, due to the quadratic dependence on the size of the larger 
dipole [See Eq. (TTT1) ]. the contribution of large values of T\ and r 2 is quite significant in the total cross section. It 
implies that in order to keep our calculations in the perturbative regime a cut in the integration on the pair separation 
should be assumed. As in [27] . we stop the r\ and r 2 integrations at a maximum dipole size, which is chosen to be 
Tmax = 1/A, with A being a free parameter. As shown in (2?|], this model successfully describes the current data 
on the total 77 cross section, on the photon structure function F 2 (x,Q 2 ) at low x and on the 7*7* cross section 
extracted from LEP doubled tagged events, with the expected value A ps Aqcd- However, that analysis did not take 
into account fluctuations effects, which is the aim of this paper. 

In what follows we will assume the IKT model for the dipole-dipole cross section and consider that the event-by- 
event scattering amplitude Af(ri,Yi) is given by Eq. ©. When taking into account the gluon number fluctuations the 
scattering amplitude will be replaced by the averaged (physical) amplitude, (J\f(p, p s )), which is given by averaging 
over all possible gluon realizations/events, corresponding to different events in an experiment, Eq. ([7]). 
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FIG. 4: (a) Energy dependence of the real <r 77 cross section, (b) Rapidity dependence of the virtual er 7 * 7 * cross section at 
Qi = Ql = Q 2 = 20 GeV 2 . 



IV. RESULTS 



We use the same values for the parameters A, xq and D obtained in Ref. [18| by fitting the F2 HERA data. 
Moreover, we assume three flavours with equal masses (my = 0.14 GeV). The free parameter in our calculations is A, 
which determines the normalization of the cross sections. As in [27| . we choose A in such a way that the experimental 
data of the real 77 cross section [47| in the low energy regime (W < 60 GeV) are well fitted. In Fig. Eta) we 
present our results for the energy dependence of the real cross section considering the IKT model without and with 
fluctuations for different values of A. We can see that it is not possible to describe the data by using the same value 
for A in the two different analyses. The values of A that allows the description of the data are A = 0.26 GeV (solid 
line) and A = 0.22 GeV (dashed line), which are near Aqcd, in agreement with our expectations. This result can 
be interpreted as an indication that the IKT model for the dipole-dipole cross section captures the main features of 
the interaction. Furthermore, we observe that the inclusion of the gluon fluctuations effects implies a smother energy 
behaviour, which agree with the theoretical expectation. However, the difference between the predictions is smaller 
than 3 % at W = 1000 GeV, which implies that tr 77 should not be the ideal observable to determine the presence of 
the gluon fluctuations effects. 

In Fig. IDJb) we present our predictions for the two-photon cross section, as a function of the rapidity Y = 
ln(W 2 /Q1Q2), for the case Q\ = Q\ (with large Q\ 2 ) corresponding to the interaction of two (highly) virtual 
photons. In this case, the contribution of the saturation effects is expected to be smaller and, consequently, gluon 
fluctuation effects to be larger. One can see that gluon number fluctuations diminish the cross section by almost 30% 
at Y = 12, and their effects keep increasing at larger values of Y. At smaller values of Q 2 , we verify that both models 
describe the experimental data with a smaller difference between them. The experimental point is taken from the L3 
Collaboration 48]. 

Fig. [5] shows our predictions for the x = Q 2 /(Q 2 + W 2 ) dependence of the photon structure function F^ix^Q 2 ) 
for different values of the photon virtualities. The basic idea is that the quasi-real photon structure may be probed 
by other photon with a large momentum transfer. We present in the lower right panel our predictions for the virtual 
photon structure function. Although there exist only very few data on this observable, its experimental study is 
feasible in future linear colliders. The current experimental data [ig], [5(| are described quite well by both models, 
with the difference between them increasing at small- a;. In particular, the difference can be of the order of 30% in the 
kinematical range which could be probed in the future linear colliders. 

In the above analysis we have used IKT model for the dipole-dipole cross section, but the description of this quantity 
is still an open question. Timneanu, Kwiecinski and Motyka [26| have proposed a different model (TKM model) to 
describe a dd in which a dd (r 1 ,r 2l Y) = a$ d N(r 2 ff ,Y). Here a$ d = (2/3)cr , with <j fixed in the GBW model by 
fitting the ep HERA data and r 2 ^ = rfr 2 /(rf + r 2 ) being an effective dipole size. The light quark mass m in the 
photon wave functions is assumed to be a free parameter, to be fixed in such a way to describe the experimental 
data for a 11 at small values of the center-of-mass energy. Although largely used in the literature, TKM model seems 
not to be well justified, for it assumes the impact-parameter factorization of the dipole-dipole cross-section, which 
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implies a oc <j$. In the dipole-proton case, ao reflects the size of the proton. However, in the dipole-dipole case, it 
should reflect the size of the larger dipole. Therefore, taking a^ d as a constant is an unphysical procedure. Thus, we 
believe that the IKT model is more realistic to describe observables in two-photon interactions. However, this subject 
deserves more detailed studies. For completeness of the present study, we present in Fig. [6] a comparison between the 
predictions without and with fluctuations obtained using the TKM model for a dd with the experimental for the 77 
cross section. As already observed before, the inclusion of the fluctuation effects implies a smother behaviour with 
energy. It is important to emphasize that the impact of the fluctuations is larger than observed in the IKT model. 
We have checked that this conclusion also is valid for the other observables. 
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FIG. 6: Comparison between IKT and TKM models in the description of the cr 77 data. 



V. CONCLUSIONS 



A current open question in the QCD dynamics at high energies is if the gluon fluctuations effects should be considered 
in the description of the observables. Results obtained using toy models indicate that these effects are suppressed by 
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the running coupling corrections to the evolution but, since the investigation in QCD of the consequences of these 
effects simultaneously is still prohibitive, phenomenological studies which test the implications of the gluon fluctuation 
effects on obscrvablcs of different processes remain important. Following previous studies that indicate that the gluon 
fluctuation effects may be present in ep collisions at HERA, in this paper we investigated their influence on some of the 
observables which could be measured in the future linear e + e~ colliders. In particular, we studied the consequences 
of fluctuations on -yM^*) interactions, in the fixed coupling case, within the dipole picture, using a dipole model for 
the dipole-dipole cross section. Our results indicate that these effects diminish the increasing with the energy of the 
real and virtual cross sections and modify the x-dependence of the photon structure function. For the total virtual 
cross section and the photon structure function, the reduction with respect to the case without fluctuations can be 
of the order of 30%, which implies that these effects should not be disregarded in the description of these observables 
in future colliders. However, because there are large uncertainties present in the current analyses — in particular, the 
normalization of cross sections, associated with the choices of quark masses and the parameter A, and even in the 
models for the dipole-dipole cross section add — observing the presence of the fluctuation effects on inclusive processes 
in photon-photon collisions will be a hard task. 
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